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Spectroscopic Evidence of Charge Exchange X-ray Emission from 
Galaxies 
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What are the origins of the soft X-ray line emission from non-AGN galaxies? XMM-Newton RGS spectra of nearby 
non-AGN galaxies (including starforming ones: M82, NGC 253, M51, M83, M61, NGC 4631, M94, NGC 2903, and the 
Antennae galaxies, as well as the inner bulge of M3 1) have been analyzed. In particular, the Ka triplet of O vn shows that 
the resonance line is typically weaker than the forbidden and/or inter-combination lines. This suggests that a substantial 
fraction of the emission may not arise directly from optically thin thermal plasma, as commonly assumed, and may instead 
originate at its interface with neutral gas via charge exchange. This latter origin naturally explains the observed spatial 
correlation of the emission with various tracers of cool gas in some of the galaxies. However, alternative scenarios, such 
as the resonance scattering by the plasma and the relic photo-ionization by AGNs in the recent past, cannot be ruled out, 
at least in some cases, and are being examined. Such X-ray spectroscopic studies are important to the understanding of the 
relationship of the emission to various high-energy feedback processes in galaxies. 
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1 Introduction 

Soft X-ray line emission has commonly been used to trace 
various types of galactic feedback in nearby starburst and 
normal galaxies. Assuming an origin of this emission in op- 
tically thin thermal (collisionally excited) hot plasma, one 
may further estimate its mass, energy, and chemical con- 
tents and even its outflow rate from such a galaxy. How- 
ever, the X-ray emission can be seriously contaminated, if 
not dominated, by various other contributions. They must be 
quantified before we can reliably use the emission to study 
the galactic feedback and its impact on the galactic ecosys- 
tem. 

Here we focus on the evidence of a significant contribu- 
tion from charge exchange (CX, also called charge transfer) 
X-ray emission at interfaces between hot and cool gases in 
galaxies. In the simplest form, a CX can be represented by 

A q+ +N^ A( q -^+* 



(1) 



where a highly ionized ion A q+ (like O vm or Ne x) picks 
up an electron from a neutral species N (like H, H2, or He), 
producing an excited ion A( q ~ x ' + *, which then decays to 
the ground state, emitting X-ray photon(s). Extensive stud- 
ies have been carried out on the solar wind CX (see Dennerl 
2010 for a recent review). In particular, the solar wind CX 
origin of the X-ray emission from comets has been estab- 
lished spectroscopically. The solar wind CX with the local 
interstellar medium may also account for a substantial part 
of the observed soft X-ray background (e.g., Koutroumpa 
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et al. 201 1). Furthermore, the CX has been speculated to be 
important in many other types of astrophysical objects (e.g., 
Lallement 2004) and has specifically been proposed to be 
responsible for emission features unaccounted for by ther- 
mal plasma modeling alone in X-ray CCD spectra of Galac- 
tic supernova remnants (e.g., Katsuda et al. 201 1) and giant 
H 11 regions (e.g., Townsley et al. 201 1). Such CX scenarios, 
if confirmed spectroscopically, could provide new insights 
into the physical, chemical, and kinematic properties of the 
underlying hot plasma and its interplay with cool gas. 

The most straight-forward X-ray spectroscopic diagnos- 
tic of a significant CX contribution is the so-called G ratio 
[= (/ + i)/r, forbidden + inter-combination to resonance 
lines] of the Ka triplets of He-like ions. Take O vn as an 
example. For an optically thin thermal plasma with a tem- 
perature of a few 10 6 K and in collisional ionization equilib- 
rium, G should be considerably less than one (Porquet et al. 
2001). The temperature may be constrained independently, 
for example, by the resonance transition ratios of H-like to 
He-like ions (e.g., Liu et al. 201 1). 

The only available instruments that allow for such X- 
ray spectroscopy of moderately extended X-ray-emitting re- 
gions are the XMM-Newton Reflection Grating Spectrom- 
eters (RGSs), which have both a large effective area and a 
large dispersion power. XMM-Newton observations, with a 
spatial resolution of ~ 6" (FWHM) at ~ 1 keV, are typ- 
ically carried out with simultaneous exposures by the two 
co-aligned RGSs, together with the European Photon Imag- 
ing Cameras (X-ray CCDs). The orientation of the satellite 
determines the dispersion direction of the RGSs, which are 
slit-less spectrometers and are sensitive to photons in the 
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0.3-2 keV range. For a point-like source, the spectral reso- 
lution of the RGSs is X/SX ~ 400 (HEW) at - 20 A. For an 
extended region, one can obtain useful 1-D spatial informa- 
tion over the cross-dispersion range of 5' wide. In the dis- 
persion direction, the RGSs cover the entire field of view of 
the telescope mirrors. The extent of the region translates to 
an apparent ( 1 st-order) wavelength shift of 0. 1 4 A arcmin - 1 
with respect to a non-displaced source. Therefore, for a re- 
gion with a moderate extent (e.g., < 2'), the spectral resolu- 
tion of an integrated RGS spectrum can still be substantially 
better than that of X-ray CCD data. In fact, the wavelength 
shift can sometimes be used to extract the spatial informa- 
tion along the dispersion direction. 

2 Examples 

2.1 M31 Bulge 

Our existing study of the central bulge of M31 (e.g., Fig.[T} 
demonstrates the spectroscopic diagnostic power of RGS 
observations in probing the nature of the soft X-ray line 
emission (Liu et al. 2010). Fig. [Ip, for example, shows a 
clear intensity excess of the observed spectrum above a sim- 
ple thermal plasma model fit at ~ 22 A. This excess is ap- 
parently due to the prominence of the forbidden line (at 22. 1 
A) with respect to the resonant one (21.6 A) of the He-like 
O vii Ka triplet. Furthermore, the excess is shown to be 
present throughout much of the cross-dispersion range of 
the RGS observations and seems to correlate well with the 
distribution of the cool gas, as traced by the Ha emission 
(Fig. QJi; Liu et al. 2010). In addition, the centroid of the 
Ne ix Ka triplet is also closer to its forbidden line wave- 
length than to its resonant one (Fig.QJ). A natural explana- 
tion of these results is that the excess represents the contri- 
bution from the CX between the two gas phases. The pres- 
ence of the diffuse hot gas in the bulge is expected from the 
stellar feedback in forms of mass loss and Type la SNe of 
low-mass stars (e.g., Tang et al. 2010). There is little evi- 
dence for any recent massive star formation in the bulge or 
significant AGN activity (Li et al. 2009, 201 1). The CX and 
other heating/evaporation processes at cool/hot gas inter- 
faces may play an important role in launching a bulge-wide, 
mass-load outflow and in starving the supermassive black 
hole at the center of the galaxy (Li & Wang 2007; Tang et 
al. 2009). Similar processes may also happen in other galac- 
tic spheroids. 

2.2 Active Star-forming Galaxies 

RGS spectra of star-forming galaxies, as well as individ- 
ual starburst regions, show roughly similar X-ray spectro- 
scopic characteristics. Fig. |2]presents two example spectra. 
Our existing analysis has focused on the O vii Ka triplet 
of a sample of nine nearby star-forming galaxies with avail- 
able XMM-Newton observations of relatively good signal to 
noise ratios (Liu et al. 2012). The observed G ratios of these 
galaxies are all larger than 1.2, substantially higher than 



expected from the thermal emission. For NGC 253, M51, 
M83, M61, NGC 4631 and the Antennae galaxies, the G ra- 
tios are consistent with the laboratory-measured ratio (2.2) 
of the CX, while for M94 and NGC 2903, a combination 
of thermal and CX gives better fits. The G ratios for other 
He-like ions such as Ne ix (with higher ionization energy 
thresholds), are typically much smaller and are thus more 
consistent with the thermal origin. In the central region of 
M82 (Fig.[3j Liu et al. 201 1), for example, the CX contribu- 
tions are 90, 50 and 30% to the fluxes of the O vii, Ne ix, and 
Mg xi Ka triplets, respectively; the averaged CX contribu- 
tion is thus ~ 50% of the total observed triplet fluxes. The 
O vii triplet G ratio also appears to change from one region 
to another. Furthermore, the cross-dispersion profiles of the 
O vii triplet and O vm Lya line are similar, indicating that 
they may have a similar origin (Liu et al. 2010). 

3 Implications 

Determining the origin(s) of the soft X-ray line emission 
is important for understanding the galactic feedback, hence 
the evolution of galaxies. Take nuclear starbursts as an ex- 
ample, which can drive so-called superwinds and have sub- 
stantial impacts on galactic ecosystem and possibly even 
on the surrounding intergalactic medium. However, it re- 
mains controversial as to whether the observed diffuse soft 
X-ray emission around starburst galaxies represents the su- 
perwinds themselves or just their interfaces with entrained 
or even in-falling cool gas clouds/filaments (e.g., Strickland 
& Stevens 2000; Bauer et al. 2007). Temperature and abun- 
dance anomalies have also been reported for starburst galax- 
ies when diffuse X-ray emission is modeled with thermal 
plasma alone. As demonstrated above, the modeling of the 
thermal emission alone is not adequate and that other con- 
tributions or effects need to be accounted for in order to 
reliably characterize the properties of the superwinds (e.g., 
Liu etal. 2011). 

If the CX is indeed significant, it will then be essential 
to quantify its contribution. In contrast to the thermal emis- 
sion, the CX contributes only lines. If an X-ray spectrum 
arising even partly from the CX is modeled with thermal 
emission only, then the inferred temperature and chemical 
properties of the plasma could be grossly wrong. To mea- 
sure the plasma properties more accurately, one needs to 
decompose the CX and thermal contributions, spectroscop- 
ically. While our existing studies are based chiefly on the 
G ratios of the Ka triplets of He-like ions (mostly O vn; 
e.g., Liu et al. 2010, 2011, 2012), jointly fitting multiple 
emission lines or even an entire spectrum will maximize the 
use of the spectroscopic information available in the RGS 
data to effectively test the CX model, to tightly constrain 
multiple model parameters with proper error propagation, 
and to decompose the thermal and CX contributions. In a 
joint analysis, even spectral lines that are not individually 
detected can be useful, as demonstrated in existing X-ray 
absoiption line spectroscopy (e.g., Yao & Wang 2005; Wang 
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Fig. 1 (a) Tri-color image of the spheroid region of M31: SpitzerfMIPS 24 fim emission (red), 2MASS K-band emission (green), and Chandra 0.5-2 
keV emission of truly diffuse hot gas (blue). The data like these show clear evidence for a subsonic outflow of hot gas (e.g., Li & Wang 2007; Li et al. 
1011; Wang 2010 and references therein), (b) Intensity contours of the diffuse soft X-ray emission overlaid on an H a image of the M31 central region. 
The contours are at 5, 10, 19, 27, 35, 45 and 55 X 10 — 3 cts s — 1 arcmin - 2 . The plus sign marks the M31. A detailed analysis of the data shows evidence 
for strong interaction between the hot and cool gas phases, which may provide a mechanism for starving the nucleus as marked by the plus sign (Li et 
al. 2009). (c) Fluxed XMM-JVewton RGS spectrum of the M31 central region (Liu et al. 2010). Both RGS1 and RGS2 spectra are combined together. 
The curves show the fit with a single-temperature thermal plasma, plus a power-law that characterizes the discrete source contribution. Notice that the 
centroid of the O VII line is closer to / wavelength than the r wavelength. [Figures are reproduced from Li et al. (2009) for (a) and (b) and from Liu et al. 
(2010) for (c).] 



et al. 2010 and references therein). Such an analysis will be- 
come possible with the upcoming CX spectral model (see 
the presentation by Smith et al. in this workshop). A joint 
analysis can, in principle, provide us with unique diagnos- 
tics of the hot plasma and its interplay with cool gas. The 
local emissivity of the CX (assuming single electron cap- 
ture) is ~ n^ q +n neu v r a, where n^ q + and n neu denote the 
ion and neutral atom densities, while v r and a are their rela- 
tive velocity and CX cross-section. The n Aq + values of rel- 
evant ions (e.g., O vm and Ne x) depend on the total parti- 
cle density, metal abundance, and ionization temperature of 
the hot plasma. These parameters can thus be inferred from 
the emission line intensities of various ion species. More- 
over, the relative intensities of multiple transitions of same 
ion species are also sensitive to v r (e.g., Greenwood et al. 
2004). Therefore, the modeling of the CX contribution can 
be a powerful tool in probing the thermal, chemical, and 
kinematical properties of the hot plasma. 

4 Potential Complications 

While the CX appears to be an important mechanism in gen- 
erating soft X-ray line emission, we also need to consider 
other plausible scenarios in various circumstances. One of 
such scenarios is the resonance scattering of emission line 
photons by the X-ray-emitting plasma itself (e.g., Porquet et 
al. 2001). Such scattering re-distributes optically thick line 
photons to outer regions. This effect on the spatial distribu- 
tion of the Fe xvn 2p — 3d transition has clearly been shown 
with the RGS observations of a few giant elliptical galaxies 
(Xu et al. 2002; Werner et al. 2009). However, the effective- 
ness of the scattering is not clear in low Lx / Lk spheroids 
(e.g., the bulge of M31), which contain hot plasma of rel- 



atively low column density, or in starburst galaxies, where 
hot superwinds are expected to have large velocity disper- 
sions. Our preliminary estimates indicate that the scatter- 
ing of the O vii Ka and O vm Lya resonance lines could 
be significant in the inner bulge region of M31, depending 
largely on how small the bulk and turbulent velocities of the 
plasma are. We are currently conducting detailed simula- 
tions to determine the significance of the resonance scatter- 
ing, which may incidentally place constraints on the plasma 
velocities. If the resonance scattering is important, then one 
may expect to observe a systematic decreasing G ratio with 
increase distances from galaxy central regions. But this ef- 
fect is not apparent in the existing data (Liu et al. 2012), 
although a thorough analysis is yet to be performed. The ef- 
fect could also be diminished if scattered photons are largely 
absorbed by cool gas. This complication also needs to be in- 
vestigated. 

Large G ratios can further be produced by plasma under 
a non-equilibrium ionization condition. Such a condition 
may be created in a cooling superwind due to fast adiabatic 
expansion, for example (e.g., Li et al. 2006). The delayed 
re-combination of highly ionized ions in the plasma tends 
to occur primarily in outer regions of a starburst galaxy. But 
the large G ratio is observed in inner regions of galaxies (Liu 
et al. 2010, 2012). The recombination scenario also does 
not explain the good correlation of the X-ray emission with 
the cool gas tracers. A similar condition, however, may also 
be produced by AGNs. Our sample selection of non-AGN 
galaxies should have minimized this complication (Liu et al. 
2012). We have even avoided the inclusion of galaxies with 
significant Seyfert 2 nuclei, which could contribute to dif- 
fuse soft X-ray emission via photo-ionization and/or scatter- 
ing. However, there might be significant contributions from 
AGN relics in some galaxies. We are modeling the emis- 
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Fig. 2 Fluxed RGS spectra of the central region of M51 (a; with a total effective exposure of 97 ks) and the Antennae galaxies (b; 172 ks). Spectral 
bins with no exposure (e.g., due to the RGS CCD gaps) are marked with bars at the wavelength axis. Positions of key lines or recombination edges are 
marked. 
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Fig. 3 (a) Fluxed RGS spectrum of M82 extracted from a central 1 ' region (Liu et al. 201 1 ; see also Ranalli et al. 2008). The continuum is mainly due 
to X-ray binaries, (b) Close-up of the O VII K« triplet (Liu et al. 2011). A continuum, which is interpolated between the line wings, has been subtracted. 
[Figures are reproduced from Liu et al. (201 1)] 



sion from diffuse plasma heated in the past, particularly as this scenario could be spectroscopically distinguished from 
a function of time since the last AGN activity to check how the CX. One known distinction is that a significant AGN 
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Fig. 4 O VIII and Ne X Lya line intensity images of M82: (a) Ne X contours overlaid on the O VIII image; (b) O VIII intensity contours overlaid on 
the Spitzer IRAC 8 fim image. The Lya line intensity distributions are reconstructed from five RGS observations. The local backgrounds (including the 
on-line continuum emissions chiefly from point-like sources) have been subtracted. 



relic should show the so-called recombination continuum 
(or edges; e.g., Guainazzi & Bianchi 2007), whereas the CX 
should not. 

In addition to the spectral analysis and modeling, spa- 
tially resolved information also needs to be fully explored. 
For example, Fig. |4] shows a preliminary version of 2-D X- 
ray line emission images of M82, reconstructed from the 
RGS observations. The comparison of such spatial distribu- 
tions of various emission lines will provide new tests of the 
scenarios. 

Ultimately, we probably also need as much information 
as can be extracted from data in other multi-wavelength ob- 
servations to understand the complex high-energy phenom- 
ena and processes involved in the heating, transportation, 
and cooling of the hot plasma and in its interaction with 
cool gas in galaxies. 
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